arXiv:1504.06030v2 [quant-ph] 22Jul2015 


Quantum theory of a bandpass Purcell filter for qubit readout 

Eyob A. Sete^*, John M. Martinis^’^, and Alexander N. Korotkov^ 
^Department of Electrical and Computer Engineering, 

University of California, 

Riverside, California 92521, USA 
^Department of Physics, University of California, 

Santa Barbara, California 93106, USA 
^ Coogle Inc., Santa Barbara, USA 
(Dated: July 23, 2015) 

The measurement fidelity of superconducting transmon and Xmon qubits is partially limited by 
the qubit energy relaxation through the resonator into the transmission line, which is also known as 
the Purcell effect. One way to suppress this energy relaxation is to employ a filter which impedes 
microwave propagation at the qubit frequency. We present semiclassical and quantum analyses 
for the bandpass Purcell filter realized by E. Jeffrey et al. [Phys. Rev. Lett. 112, 190504 (2014)]. 
For typical experimental parameters, the bandpass filter suppresses the qubit relaxation rate by up 
to two orders of magnitude while maintaining the same measurement rate. We also show that in 
the presence of a microwave drive the qubit relaxation rate further decreases with increasing drive 
strength. 
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I. INTRODUCTION 

The implementation of fault-tolerant quantum in¬ 
formation processing [I| requires high-fidelity quantum 
gates and also needs sufficiently fast and accurate qubit 
measurement. Superconducting quantum computing 
technology inn is currently approaching the threshold 
for quantum error correction. Compared with the recent 
rapid progress in the increase of single-qubit and two- 
qubit gate fidelities, qubit measurement shows somewhat 
slower progress. The development of faster and higher- 
fidelity qubit readout remains an important task. 

In circuit QED (cQED) [lllllll) the qubit state is in¬ 
ferred by measuring the state-dependent frequency shift 
of the resonator via homodyne detection. This method 
introduces an unwanted decay channel [l^ for the qubit 
due to the energy leakage through the resonator into the 
transmission line, the process known as the Purcell effect 
0[il. The Purcell effect is one of the limiting factors 
for high fidelity qubit readout. 

In principle, the Purcell rate can be suppressed by 
increasing the qubit-resonator detuning, decreasing the 
qubit-resonator coupling, or decreasing the resonator 
bandwidth due to damping. However, these simple meth¬ 
ods increase the time needed to measure the qubit. This 
leads to a trade-off between the qubit relaxation and mea¬ 
surement time, whereas it is desirable to suppress the 
Purcell rate without compromising qubit measurement. 
Several proposals have been put forward for this pur¬ 
pose, which include employing a Purcell filter |16l - ll^ . 
engineering a Purcell-protected qubit [l^, HH, or using 


‘Present address: Rigetti Quantum Computing, 2855 Telegraph 
Ave, Berkeley, CA 94705, USA 


a tunable coupler that decouples the transmission line 
from the resonator during the qubit-resonator interac¬ 
tion, thereby avoiding the Purcell effect altogether . 

The general idea of the Purcell filter is to impede the 
propagation of the photon emitted at the qubit frequency, 
compared with propagation of the microwave field at the 
resonator frequency, used for the qubit measurement. A 
notch (band-rejection) filter detuned by 1.7 GHz from 
the resonator frequency was realized in Ref. [l^. A fac¬ 
tor of 50 reduction in the Purcell rate was demonstrated 
when the qubit frequency was placed in the rejection 
band of the filter. A bandpass filter with the quality 
factor Q{ ~ 30 (and corresponding bandwidth of 0.22 
GHz) centered near the resonator frequency was used in 
Ref. 0. This allowed the qubit measurement within 
140 ns with fidelities Ep) = 98.7 and F|o) = 99.3 for the 
two qubit states. (The bandpass Purcell filter was also 
used in Ref. [l^ ; it had a similar design with a few minor 
changes.) A major advantage of the bandpass Purcell fil¬ 
ter in comparison with the notch filter is the possibility to 
keep strongly reduced Purcell rate for qubits with prac¬ 
tically any frequency (except near the filter frequency), 
thus allowing quantum gates based on tuning the qubit 
frequency, and also allowing multiplexed readout of sev¬ 
eral qubits by placing readout resonators with different 
frequencies within the filter bandwidth. 

In this work, we analyze the Purcell filter of Ref. [13 
using both semiclassical and quantum approaches and 
considering both the weak and the strong drive regimes. 
Our semiclassical analysis uses somewhat different lan¬ 
guage compared to the analysis in Ref. [l^]; however, 
the results are very similar (they show that with the 
filter the Purcell rate can be suppressed by two orders 
of magnitude, while maintaining the same measurement 
time). The results of the quantum analysis in the regime 
of a weak measurement drive or no drive (considering 
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the single-photon subspace) practically coincide with the 
semiclassical results. In the presence of strong microwave 
drive, the Purcell rate is further suppressed with increas¬ 
ing drive strength. We have found that this suppression 
is stronger than that obtained without a filter (23l |. 

In Sec. II we discuss the general idea of the bandpass 
Purcell filter and analyze its operation semiclassically. 
Section III is devoted to the quantum calculation of the 
Purcell rate in the presence of the bandpass Purcell filter. 
In Sec. IV we discuss further suppression of the Purcell 
rate due to an applied microwave drive. Section V is the 
Conclusion. In the appendix we review the basic theory 
of a transmon/Xmon qubit measurement, the Purcell de¬ 
cay, and the corresponding measurement error without 
the Purcell filter. 


II. IDEA OF A BANDPASS PURCELL FILTER 
AND SEMICLASSICAL ANALYSIS 


In the standard cQED setup of dispersive measurement 
(Fig. HI) the qubit interaction with the resonator slightly 
changes the effective resonator frequency depending on 
the qubit state, so that it is when the qubit is in the 
excited state and when the qubit is in the ground 
state. The dispersive coupling x is defined as 

( 1 ) 


In the two-level approximation for the qubit, y = 
— wj’), where g is the qubit-resonator coupling 
and Wq and wj? are the bare frequencies of the qubit 
and the resonator, respectively |ll| . For a transmon 
or an Xmon qubit, y is usually significantly smaller, 
y Ri - w)")], where 5q is the 

qubit anharmonicity (Jq > 0); moreover, y as well as the 
central frequency {iJf' -\- a;i®^)/2 depend on the number 
of photons n in the resonator (see and the Ap¬ 

pendix for a more detailed discussion). The resonator 
frequency change (and thus the qubit state) is sensed 
by applying the microwave field with amplitude e, then 
amplifying the transmitted or reflected signal, and then 
mixing it with the applied microwave field to measure its 
phase and amplitude (Fig. [T|). 

In the process of measurement, the qubit decays with 
the Purcell rate m 


r 


(Wq - UJrY 


( 2 ) 


where k is the resonator energy damping rate (mostly 
due to leakage into the transmission line - see Fig. [IJ. 
Note that in this formula we do not distinguish the bare 
and effective frequencies. In the quantum language this 
can be interpreted as the leakage with the rate k of the 
qubit “tail” g'^/{ujq — Wj-)^, existing in the form of the 
resonator photon. However, the Purcell decay also has a 
simple classical interpretation via the resistive damping 
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FIG. 1: Schematic of a standard circuit QED qubit readout 
setup. The qubit state slightly changes the resonator fre¬ 
quency LJi (due to qubit-resonator interaction with strength 
g), and this is sensed by passing the microwave through (or 
reflecting from) the resonator. The amplified outgoing mi¬ 
crowave is combined with the local oscillator at the mixer, 
whose output is measured to discriminate the qubit states. 
The energy decay k of the resonator is mainly due to its cou¬ 
pling with the transmission line. 


[T^, essentially being a linear effect, in contrast to the 
dispersion (HD - see Appendix for more details, including 
dependence of P on n [23| . 

The Purcell decay leads to measurement error; there¬ 
fore, it is important to reduce the rate P. This can be 
done by decreasing the ratio g/\ujq — WrI; however, this 
decreases y and thus increases the necessary measure¬ 
ment time tm (see Appendix for more details). Another 
way to decrease P is to use a very small leakage rate k; 
however, this also increases the measurement time tm be¬ 
cause the ring-up and ring-down processes give a natural 
limitation tm ^ and in many practical cases it is 

even t„i 10 

It would be good if the rate k which governs the mea¬ 
surement time were different from k in Fq. specifi¬ 
cally if K for the Purcell decay were much smaller than 
K for the measurement. This is exactly what is achieved 
by using the bandpass filter of Ref. [l^ ■ There are other 
ways to explain how this Purcell filter works [13 , but 
here we interpret the main idea of the bandpass Pur¬ 
cell filter as producing different effective rates Keff for 
the measurement and for the Purcell decay (so that the 
measurement microwave easily passes through the filter, 
while the propagation of the photon emitted by the qubit 
is strongly impeded by the filter). 

The schematic of the qubit measurement with the 
bandpass Purcell filter of Ref. [13 is shown in Fig. [21 Be¬ 
sides the readout resonator with qubit-state-dependent 
frequency Wr = or Wj. = (jjf\ there is a second (fil¬ 
ter) resonator with frequency Wf, coupled with the read¬ 
out resonator with the coupling Q. (The coupling Q is 
inductive, but we draw it as capacitive to keep the figure 
simple.) The filter resonator leaks the microwave into the 
transmission line with a relatively large damping rate Kf, 
so that its Q-i&cioT is Qf = Wf/^f ~ 30, while \G\ ^ Kf- 
The leaked field is then amplified and sent to the mixer 
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FIG. 2: Qubit measurement schematic with the bandpass 
Purcell Hlter of Ref. The readout resonator with fre¬ 

quency LJi (which depends on the qubit state) is coupled (cou¬ 
pling Q) with a filter resonator of frequency uj{, which decays 
into the transmission line with the rate K{. The further pro¬ 
cessing of the outgoing microwave (not shown) is the same 
as in Fig. [T] The microwave drive can be applied either to 
the readout (sr) or to the filter (sf) resonators. Coupling 
with the decaying Hlter resonator produces an effective de¬ 
cay rate Kefi of the readout resonator, which depends on the 
drive frequency. As a result, for the measurement microwave 
Keff = Kr, while the qubit sees a much smaller value Kefr = ACq, 
thus leading to a suppression of the qubit Purcell decay by a 
factor Kq/ftr. 


(not shown) in the same way as in the standard cQED 
setup. The readout and filter resonators are in general 
detuned from each other, but not much, |wi. — Wfj ^ K{ 
(detuning is needed to multiplex readout of several qubits 
using the same filter resonator 0113; for simplicity we 
consider the measurement of only one qubit). The filter 
resonator is pumped with the drive frequency Wd (close 
to Wr) and amplitude Sf. However, for us it will be eas¬ 
ier to first assume instead that the readout resonator is 
pumped with amplitude Sr (Fig. [2|), and then show the 
correspondence between the drives and Sf. 

Let us use the rotating wave approximation [2^, 
with the rotating frame based on the drive fre¬ 

quency Wd. Then the evolution of the classical field am¬ 
plitudes a{t) and f3{t) in the readout and filter resonators, 
respectively, is given by the equations 

a = —iArdCt — iQP — (3) 

/3 = -*Afd/3-*e*a-y/3, (4) 

where a and (3 are normalized so that jap and |/3p are 
the average number of photons in the resonators, Sr is 
normalized correspondingly, and 

Ard = Wr - Wd, Afd = Wf - Wd (5) 


(recall that Wr depends on the qubit state). If we are 
not interested in the details of evolution on the fast time 
scale then we can use the quasisteady state for /3 
[obtained from Eq. (jH) using /3 = 0], 


/3 


-ig* 

- TTi -1- Cl, 

Ki/2 -\-1 Afd 


( 6 ) 


which can then be inserted into Eq. ®, giving 
d = —i(Ard + Sujr) a -^ a — iSr, 

4\g\^ 1 

Atf l+(2Afd/Atf)2’ 

^ _ |t/pAfd Afd 

iKi/2r + Aj,-- Ki 

Thus we see that the field a evolves in practically the 
same way as in the standard setup of Fig. |T] however, 
interaction with the filter resonator shifts the readout 
resonator frequency by Sujr and introduces the effective 
leakage rate KeS of the readout resonator. 

Most importantly, Keff depends on the drive frequency. 
For measurement we use Wd ~ Wr, so KeS is approximately 


( 7 ) 

( 8 ) 
(9) 




4|gP 

K{ 


1 

1 + [2(a;r - ujf)/KiY 


( 10 ) 


However, when the qubit tries to leak its excitation 
through the readout resonator, this can be considered 
as a drive at the qubit frequency, Wd = Wq, and the cor¬ 
responding Atefi is then 


_ 4|gp _1_ 

“ Kf 1 -F [2(a;q - uj{)/KiY ’ 


( 11 ) 


which is much smaller than Kr if the qubit is detuned 
away from the filter linewidth, jwq — Wfj Kf. This 
difference is exactly what we wished for suppressing the 
Purcell rate T: the measurement is governed by Kr, while 
the qubit sees a much smaller value Kq. Therefore, we 
would expect that the Purcell rate is given by Eq. @ 
with K = Kq [see Eq. (1321) later], while the “separation” 
measurement error is given by Eqs. (IA36|) ~ (|A38|) with 
K = Kr (see Appendix). As a result, compared with the 
standard setup (Fig. [1]) with the same physical parame¬ 
ters for measurement, the Purcell rate is suppressed by 
the factor 

Kq 1+[2{UJ,-Uj{)/Kf]^ 

Kr l+[2(cdq-Cdf)/Kf]2 ^ ^ 


This is essentially the main result of this paper, which 
will be confirmed by the quantum analysis in the next 
section. (To avoid a possible confusion, we note that Kq 
is not the qubit decay rate; it is the resonator decay, as 
seen by the qubit.) 

Our result for the Purcell suppression factor was based 
on the behavior of the field amplitude in the readout 
resonator. Let us also check that the field yti propagating 
in the outgoing transmission line behaves according to 
the effective model as well. The outgoing field amplitude 
is 7ti = \/Kf/3 (in the normalization for which jytip is 
the average number of propagating photons per second). 
Using Eq. (|6]), we find 


7ti 


Kjf j2 -)- zAfrj 


a, 


(13) 
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so, as expected, the outgoing amplitude behaves as in 
the standard setup of Fig. [T] with k — Keffj up to an 
unimportant phase shift (p = arg[—iC/*/(Kf/2 + iAfd)]. 
Note that to show the equivalence between the dynamics 
(including transients) of the systems in Figs. [T] and [2] 
we needed the assumption of a sufficiently large Kf in 
order to use the quasisteady state (0). However, this 
assumption is not needed if we consider only the steady 
state (without transients). 

So far we assumed that the measurement is performed 
by driving the readout resonator with the amplitude 
Now let us consider the realistic case [iS[i3 when the 
drive Ef is applied to the filter resonator. The evolution 
equations m and (U) for the classical field amplitudes 
are then replaced by 

a = —zArdO — iGP, (14) 

/3 =-iAfdlS - ig*a - Y l3 - ( 15 ) 

so that the quasisteady state for the filter resonator is 


the transfer function when driving the readout resonator, 
the transfer function for driving the filter is (in the steady 
state) 


JtT ^ \/^ _ 2Ard/KeS 

£f Kil2 + zAfd ^ 2z(Ard + ^^r) ’ 

^eff 

where Keff can be replaced with k^. (Note a non-standard 
normalization of the transfer function because of differ- 

/r\ 

ent normalizations of and £f.) This line shape for 
the amplitude |7)j /£f| shows a dip near Wr (note that 
= 0 at Wd = Wr) and is significantly asymmetric 
when Soj-c is comparable to up, this occurs when the de¬ 
tuning between the readout and filter resonators is com¬ 
parable to Kf - see Eq. In terms of the held a in the 
readout resonator, the outgoing held at steady state is 


(f) ^ Aj-d 

5'ti =-^- 


( 20 ) 


/3 = , a + , , (16) 

Kf/2-|-zAfd Kf/2-|-zAfd 

and the held evolution in the readout resonator is 

o = -.(A,J + fo.) O - ^ a - I eu (17) 

with the same KeB and given by Eqs. ([8]) and (|9]). The 
only difference between the effective evolution equations 
© and (113 is a linear relation, 

£r o-zeft//(Kf/2-I-zAfd), (18) 

between the drive amplitudes £r and £f producing the 
same effect. Therefore, our results obtained above remain 
unchanged for driving the hlter resonator, and the Purcell 
rate suppression factor is still given by Eq. ( 113 - 
Note that in the quasisteady state the separation be¬ 
tween the hlter amplitudes /3 for the two qubit states 
does not depend on whether the drive is applied to the 
hlter or readout resonator, as long as we use the corre¬ 
spondence da between the drive amplitudes. The same 
is true for the separation between the outgoing helds yti. 
Similarly, the separation between the outgoing helds for 
the two qubit states is the same (up to the phase (p) as 
in the standard setup of Eig. [T] with £ = £r, k = k^, 
and the resonator frequency adjusted by 6uJr given by 
Eq. ©. Therefore, these conhgurations are equivalent 
to each other from the point of view of quantum mea¬ 
surement, including interaction between the qubit and 
readout resonator, extraction of quantum information, 
back-action, etc. 

Nevertheless, driving the hlter resonator produces a 
different outgoing held yti = ydtf/3, which now contains 
an additional term —i£{y/^/{K{/2 + iA{d) in comparison 
with Eq. m, which comes from the second term in Eq. 
(1161) . In particular, instead of the Lorentzian line shape of 


instead of Eq. (HI for driving the readout resonator. 

The difference between the outgoing helds and 
when driving the hlter or readout resonator (for the same 
a, i.e., the same measurement conditions) may be impor¬ 
tant for saturation of the microwave ampliher. The ratio 
of the corresponding outgoing powers is 


l7n¥ ^ f 4 

\ Kr / 1-1-(2Afd/Kf)2’ 


( 21 ) 


where we assumed |Ard| -C Kf (so that Keff ~ Kr). 
For example, if the drive frequency is chosen as Wd = 
(u!p^ + Ljpp/2, then A^d = ±x; if in this case |x| <C Kj., 
then driving the hlter resonator is advantageous because 
it produces less power to be amplihed, while driving the 
readout resonator is advantageous if |x| ^ Kj.. However, 
when Wd ^ (wi®^ -I-a;i'^^)/2, the situation is more compli¬ 
cated. 

Figure[3]shows the transient phase-space evolution /3(t) 
of the held (coherent state) in the hlter resonator when 
a step-function drive is applied to the readout resonator 
(red curves) or to the hlter resonator (blue curves), with 
the drive amplitudes related via Eq. m, and for real 
£r. The solid curves show the evolution when the qubit 
is in the excited state, while the dashed curves are for 
the qubit in the ground state. We have used parame¬ 
ters similar to the experimental parameters of Ref. [13: 

/2tx = 6.8 GHz, 127: = 6.803 GHz (so that 
2x/27r = -3 MHz), a;f/27r = 6.75 GHz, Qf = 30 (so that 
= 0.71 ns), and k“^ = 30 ns (so that Q/2Tr = 18.9 
MHz). The held evolution is calculated using either Eqs. 
©"O or Eqs. (O-dTSl); in simulations we have ne¬ 
glected the dependence of x and (wi^^ -I- )/2 on the 

average number of photons hr = \ap in the readout res¬ 
onator (see Appendix). The black dots indicate the time 
moments every 10 ns between 0 and 100 ns, and then 
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FIG. 3: Phase-space transient evolution of qubit-state- 
dependent coherent states in the filter resonator for driving 
either the readout resonator (red curves, left) or the filter 
(blue curves, right), with the drive amplitudes related via Eq. 
(na. Solid curves are for the qubit state |e), dashed curves 
are for the state \g). See text for the assumed parameters. For 
panel (a) we choose the drive frequency uJd symmetrically for 
the readout resonator, so that = 50. For panel (b) 

we choose uJd symmetrically for the filter resonator, so that 
when driving the filter; we choose = 50. The 
black dots indicate time moments in the evolution every 10 ns 
until 100 ns, then every 50 ns. Circles illustrate the coherent 
state error circles in the steady state. 


every 50 ns. The circles illustrate the coherent state er¬ 
ror circles in the steady state. We see that when 
the drive is applied to the filter resonator (blue curves), 
the evolution /3(t) is initially very fast (governed by Kf), 
while after the quasisteady state is reached, the evolution 
is governed by a much slower eventually approaching 
the steady state. When the drive is applied to the readout 
resonator (red curves), the transient evolution is always 
governed by the slower decay k^- 

In Fig. Ha) we choose the drive frequency Wd symmet¬ 
rically from the point of view of the readout resonator 
field, so that in the steady state = 50, where 


and are the average photon numbers for the two 
qubit states. For that we need Wd = +wi®^)/2-|-5wr 

with SuJr given by Eq. (H; for our parameters 510-^1 2t: = 
1.23 MHz, so uJdl^TT = 6.80273 GHz. Such symmetric 
choice of the drive frequency provides the largest sepa¬ 
ration between the two coherent states for a fixed drive 
amplitude if 2 |x| < k^- While the field in the read¬ 
out resonator is always symmetric in this case. Fig. H^-) 
shows that the field in the filter resonator is symmet¬ 
ric only when driving the readout resonator (red curves, 
= 1.2), while it is strongly asymmetric when 
the filter resonator is driven (blue curves, fi[®^ = 0.01, 
= 1.0; h[®^ is very small because for our parame¬ 
ters SuJr ~ Ixl and therefore Wd ~ wi®^). The number 
of photons in the filter is much less than in the readout 
resonator because Ki ^ k^. 

In Fig. Hb) we choose Wd so that in the steady state 
h[®^ = for driving the filter; this is the natural choice 
for decreasing the microwave power to be amplified. This 
occurs at = 6.80120, which is close to the expected 

value (wi®^ +uif^)/2, but not equal because of the asym¬ 
metry of the line shape (ITOl) . We choose the amplitudes 
to produce n^f' = 50 (then fil®^ = 22). The difference 
between and h[®^ leads to different values =1.2 

and hf®^ = 0.5 when driving the readout resonator, while 
for driving the filter the field in the filter is symmetric, 
n^f' = h[®^ = 0.2. Compared to the case of Fig. Hn); 
there is 5 times less power to be amplified for the |e) state 
(when driving the filter); however, the state separation is 
1.3 times smaller (in amplitude) for the same Thus, 
there is a trade-off between the state separation and am¬ 
plified power in choosing the drive frequency. Comparing 
the red and blue curves in Fig. Hb), we see that in the 
steady state ri[®^ and are smaller for driving the fil¬ 
ter rather than the readout resonator. This is beneficial 
because there is less power to be amplified; however, the 
ratio is not very big (as expected for a moderate value 
|X|K = 0.28). 

Note that our definition of Kr in Eq. (HOD is not strictly 
well-defined because the resonator frequency Wr depends 
on the qubit state, and the drive frequency can be in be- 
|e) Is) 

tween ' and w) '. However, this frequency difference 
is much smaller than K{, and therefore not important for 
practical purposes in the definition of k^. In an experi¬ 
ment Kr can be measured either via the field decay 0 
or via the linewidth of the steady-state transfer function 
showing the dip of | 7 ^y/ef| near the resonance Wd = Wr 
[Eq. ([HD], since near the dip KeS ~ k^- 

Thus far we assumed that all decay K{ in the filter 
resonator is due to the leakage into the outgoing 
transmission line. If k™* < K{ and the decay K{ — 
is due to leakage into the line delivering the drive Sf or 
due to another dissipation channel, then the only dif¬ 
ference compared to the previous discussion is the ex¬ 
tra factor for the outgoing field yti. This will 
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lead to multiplication of the overall quantum efficiency 
of the measurement by and will only slightly af¬ 

fect the measurement fidelity. Adding dissipation in the 
readout resonator with rate Kr,d increases the effective 
linewidth to KeS -I- Kr,d and multiplies the quantum ef¬ 
ficiency by + Kr,d). Most importantly, since Kr,d 

does not change with frequency, the Purcell suppression 
factor m becomes (Kq -I- «:r,d)/(Kr -|- Kr,d), so that the 
Purcell filter performance deteriorates; we will discuss 
this in a little more detail in Sec. IIIC. 

Note that our main result (HU for the Purcell sup¬ 
pression factor is slightly different from the result F = 
[/vf/2(a;q — a;r)]^(wq/a;r)^, which was derived in Ref. 0 
using the circuit theory. The reasons are the following. 
First, in the derivation of [l^ it was assumed that the 
two resonators have the same frequency, which makes the 
numerator in Eq. (1121) equal to 1. Second, the term 1 in 
the denominator in Eq. (IT^ was essentially neglected in 
comparison with the larger second term. Finally, the role 
of the factor (wq/wr)^ is not quite clear. In the deriva¬ 
tion of Ref. keeping this factor was exceeding the 
accuracy of the derivation, while in our derivation we 
essentially use the rotating wave approximation, which 
assumes Wq/wr ~ 1. Aside from these small differences, 
our result coincides with the result of Ref. [i3- 

III. QUANTUM ANALYSIS IN 
SINGLE-EXCITATION SUBSPACE 

In this section we discuss the quantum derivation of the 
Purcell rate in the presence of the bandpass filter in the 
regime when the resonators are not driven or driven suf¬ 
ficiently weakly to neglect dependence of the Purcell rate 
on the number of photons in the resonator . More pre¬ 
cisely, we consider the quantum evolution in the single¬ 
excitation (and zero-excitation) subspace. We apply two 
methods: the wavefunction approach, in which we use 
a non-Hermitian Hamiltonian with a decaying wavefunc¬ 
tion, and the more traditional density matrix analysis. 

In the absence of the drive and in the rotating wave 
approximation, the relevant Hamiltonian of the system 
shown in Fig. [5] (without considering decay K{) is {h = 1) 

H = WqCT+CT- -I- oj^a'a + -|- g(aV_ -I- aa+) 

+ga^b + g*ab\ (22) 

where is the bare qubit frequency, a;(’ is the bare fre¬ 
quency of the readout resonator, Wf is the filter resonator 
frequency, raising/lowering operators cr+ and cr_ act on 
the qubit state, and a are the creation and annihilation 
operators for the readout resonator, fof and b are for the 
filter resonator, g is the qubit-readout resonator coupling, 
and g is the resonator-resonator coupling. For simplicity 
we assume a real positive g, but g can be complex for 
generality (for the capacitive or inductive coupling be¬ 
tween the resonators, g is real if the same generalized 
coordinates are used for both resonators). 


Note that in the case without drive it is sufficient 
to consider only two levels for the qubit because only 
the single-excitation (and zero-excitation) subspace is in¬ 
volved in the evolution, and therefore the amount of qubit 
nonlinearity due to the Josephson junction is irrelevant. 
However, in the presence of a drive (considered in the 
next section) it is formally necessary to take into account 
several levels in the qubit (as done in the Appendix). 
Nevertheless, to leading order the Purcell rate is insensi¬ 
tive to this, because the Purcell decay is essentially a clas¬ 
sical linear effect (see discussion in the Appendix). Also 
note that the lab-frame Hamiltonian dUD assumes the 
rotating wave approximation (as in the standard Jaynes- 
Cummings Hamiltonian), since it neglects the “counter¬ 
rotating” terms of the form aV+, acr_, a^b^ , and ab. This 
requires assumption that jwq — wj’j, |wf — wj?!, 5 , and \g\ 
are small compared to w)?. 

Let us choose the rotating frame with frequency Wq, 
i.e., Hq = a;q((T+cr_ -I- a^o -1- b^b); then the interaction 
Hamiltonian V = H — Hq is 

V = ArqO^a -I- Afqb^b + g{a^a- + aa+) + ga)b + g*ab\ 

(23) 

where 

A,q = a;J>-a;^, Afq = Wf - (24) 

and the interaction picture is equivalent to the 
Schrodinger picture because exp(iiLot)H exp(—ziLot) = 
V, which is because the starting Hamiltonian (1^^ al¬ 
ready assumes the rotating-wave approximation. The 
master equation for the density matrix p, which includes 
the damping K{ of the filter resonator is 

p = —i\y, p] + Kf — b^bp/2 — pb^b/2) . (25) 

In general, the bare basis is \jnm), where j represents 
the qubit states, while n and m represent the readout and 
filter resonator Fock states, respectively. However, in this 
section we consider only the single-excitation (and zero- 
excitation) subspace, so only four bare states are relevant: 
|e) = |e 00 ), |r) = |glO), |f) = |g 01 ), and |g) = |g 00 ). 

Note that the interaction hybridizes the bare states of 
the qubit and the resonators. (Hybridization of the read¬ 
out resonator mode is essentially what makes the qubit 
measurement possible.) Therefore, when discussing the 
Purcell rate for the qubit energy relaxation, we actually 
mean decay of the eigenstate, corresponding to the qubit 
excited state. This makes perfect sense experimentally, 
since manipulations of the qubit state usually occur in 
the eigenbasis (adiabatically, compared with the qubit 
detuning from the resonator). 

A. Method I: Decaying wavefunction 

Instead of using the traditional density matrix lan¬ 
guage for the description of the Purcell effect [l^, it is 
easier to use the language of wavefunctions, even in the 
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presence of the decay Kf [2^. Physically, the wave func¬ 
tions can still be used because in the single-excitation 
subspace unraveling of the Lindblad equation corre- 
sponds to only one “no relaxation” scenario (see, e.g., 
[29|), and therefore the wavefunction evolution is non¬ 
stochastic. Another, more formal way to introduce this 
language, is to rewrite the master equation (El]) as [13, HH 
p = —i[77eff,/o] + Kfbpb"^, where Hes = V — iKfb^b/2 is 
an effective non-Hermitian Hamiltonian. Next, the term 
Kfbpb^ can be neglected because in the single-excitation 
subspace it produces only an “incoming” contribution 
from higher-excitation subspaces, which are not popu¬ 
lated. Therefore, in the single-excitation subspace we 
can use p = —i[HeSj p]- Equivalently, {ip) = —iHesltp), 
which describes the evolution of the decaying wavefunc¬ 
tion \tp{t)) = Ce(t)|e) + Cr{t)\r) + Cf(t)|f). Therefore, the 
probability amplitudes Ce,r,f satisfy the following equa¬ 
tions: 

Ce = -igcr, (26) 

Cr = —iArqCr — igCe — iGCf, (27) 

Cf = -lAfqCf - iC/*Cr - (Kf/2) Cf, (28) 

while the population pgg of the zero-excitation state |g) 
evolves as pgg = Kf|cffp or can be found as pgg = 1 — 
|ceep —|crri^—|cffp. Note that Eqs. (I?7l) and (1^ exactly 
correspond to the classical equations m and 0 with 
Ard replaced with A^q, also Afd replaced with Afq, and 
Er replaced with gce- 

Erom the eigenvalues Ae,r,f = re,r.f/2 of the 

matrix representing Eqs. (I26P "(I^5 |) . one can obtain the 
eigenfrequencies a;e,r,f and the corresponding decay rates 
Ee.r.f- These eigenvalues can be found from the qubic 
equation 

-|- A^(iArq -|- iAfq -|- Kf/2) -|- A( —ArqAfq -|- \G\'^ + g^ 

-f *Ai.qKf/2) -|- ^^(zAfq -|- Kf/2) = 0. (29) 

We are interested in the Purcell rate P = Pe, which corre¬ 
sponds to the decay of the eigenstate close to |e). Since 
Ae is close to zero, in the first approximation we can 
neglect the term A^ in Eq. (uni), thus reducing it to the 
quadratic equation. If more accuracy is needed, the equa¬ 
tion can be solved iteratively, replacing A^ with the value 
found in the previous iteration (the second iteration is 
usually sufhcient). 

Besides finding the Purcell rate P exactly or approxi¬ 
mately from Eq. (1291) . we can find it approximately by 
using quasisteady solutions of Eqs. (EH) and ( 1 ^ . to 
a large extent following the classical derivation in the 
previous section. Assuming Cf = 0 in Eq. (I28p . we 
find Cf = —iQ*Cr/{iA{q + K{f2). Inserting this quasis¬ 
teady value into Eq. EB and assuming Cr = 0 , we find 
Cr = —* 5 Ce/[iAi.q-|-|C/p/(iAfq-|-Kf/2)]. Substituting this 
quasisteady value into Eq. (l25l) . we obtain 

= "zA^q-l-lelvW + WS) 


Finally, we obtain the Purcell rate as P = —2Re(Ae), 

r =__ (31) 

A?J(Afq-|C;|VArq)2 + («f/2)2] ^ ^ 

_ ^ /or,', 

^A^^[Al + {Kt/2y] A?/ ^ ^ 

where is given by Eq. m and we assumed |C/p <C 
AfqArq to transform Eq. (ISTI) into Eq. (l32l) . 

The Purcell rate given by Eq. (I32p is exactly what we 
expected from the classical analysis in Sec. II: in the usual 
formula El) we simply need to substitute k with the read¬ 
out resonator decay rate Kq seen by the qubit. Since the 
measurement is governed by a different decay rate Kr , the 
effective Purcell rate suppression factor is given by Eq. 
EB, as was expected. This confirms the results of the 
classical analysis in Sec. II. 


B. Method II: Density matrix analysis 

We can also find the Purcell rate in a more traditional 
way by writing the master equation (I25p explicitly in the 


single-excitation subspace: 

Pee '^giPer Pre) ■> 

Pev — ^AqrPer ig{Prr Pee') F iQ /^ef, (^^) 

Pef — ^Pef ^AqfPef F iQPer ipPrf^ (^^) 

Prr — ig{Per Pre) iGPfr F iG Prf: 

Prf = “Y Prf F iAfrPrf “ iG Pf£ F iG* Prr “ ipPef, (37) 
Pff=-KfPff-iG*Prf + iGpfr- (38) 


Note that pgg = Kipff and pee F Prr F Pff F Pgg = I (how¬ 
ever, we do not use these two equations for the derivation 
of the Purcell rate). 

Using the quasisteady solutions of Eqs. El-Ell), i-e. 
assuming per = Pef = Prr = Prf = Pff = 0, we can 
obtain a lengthy equation for per, which is proportional 
to Pee. If we use the first-order expansion of this equation 
in the coupling g and neglect g^ terms (there is no g^ 
contribution), then 


_ _ 11 _ /qO', 

*ArqF|e|V(*AfqF«f/2)^""’ 

which has the form similar to Eq. El- Substituting 
this Per into Eq. (I33L we obtain the evolution equation 
Pee = —Ppee with P given exactly by Eq. El- If we 
do not use the above-mentioned approximation for the 
quasisteady per, then the result for the Purcell rate is 
slightly different and much lengthier, 

P = p^iap^f [(A,qAfq - F (A?q + p^)(«f/2)2 

F5'^(Afq -I- 2AfqArq — |C/p) -I- p'*] . (40) 







Thus the derivations based on the wavefunction and den¬ 
sity matrix languages using the quasisteady-state approx¬ 
imation both lead to practically the same result for the 
Purcell rate P. The most physically transparent result is 
given by Eq. (I32]), which corresponds to the semiclassical 
analysis in Sec. II and simply replaces k in Eq. ([2|) with 
Kq seen by the qubit, in contrast to the measurement 
process, which is governed by Kr- 

As an example, let us use the parameters similar to 
that in Ref. [l^: a;q/27r = 5.9 GHz, Wr/27r = 6.8 GHz, 
uj{j2TT = 6.75 GHz, Q{ = 30 (so that = 0.71 ns), 
g/2Tr = 90 MHz, and k~^ = 30 ns (so that G/2.tt = 18.9 
MHz). In this case the resonator decay [Eq. (ITTI) ] seen by 
the qubit is = 1.45/rs, the Purcell rate [Eq. (15^ ] is 
P = 1/(145/xs^ and the Purcell rate suppression factor 
[Eq. ([HI)] is F = 30 ns/1.45/xs= (1-b 0.442)/(l7.6^) = 
0 . 021 . 

Thus, for typical parameters the bandpass Purcell filter 
suppresses the Purcell decay by a factor of ~50. It is 
easy to increase this factor to 100 by using a;q/27r = 5.5 
GHz in the above example; however, further decrease 
of the Purcell rate is not needed for practical purposes, 
while increased resonator-qubit detuning decreases the 
dispersive shift 2% (in the above example 2x/27r Ri —3 
MHz for the qubit anharmonicity of 180 MHz, while for 
LOql2'K = 5.5 GHz the dispersive shift becomes twice less). 

Note that for the parameters in the above example, Eq. 
(1321) overestimates the exact solution for P via Eq. dMl) 
by 5%, the same 5% for Eq. (I3T|) . and Eq. (1401) overesti¬ 
mates the Purcell rate by 2%. The solution of Eq. (1^ as 
a quadratic equation neglecting gives P, which over¬ 
estimates the exact solution by 22%, while the second 
iteration is practically exact (—0.01%). The inaccura¬ 
cies grow for smaller resonator-qubit detuning (crudely 
as but remain reasonably small in a sufficiently 

wide range; for example, Eq. (1321) overestimates the Pur¬ 
cell rate by 50% for a;q/27r = 6.5 GHz, i.e. detuning of 
0.3 GHz. 


C. Nonzero readout resonator damping 

In the quantum evolution model (1251) we have consid¬ 
ered only the damping of the filter resonator with the 
rate K{. If there is also an additional energy dissipation 
in the readout resonator with the rate Kr.d (e.g., due to 
coupling with the transmission line delivering the drive 
ffr in Fig. [5]), then the master equation (1^ should be 
replaced with 

p = — i[V, p] + Kf {bpb^ — l)^hp/2 — ph''h/2) 

+ Kj^^diapa' — a‘ap/2 — pa^a/2). (41) 

In the wavefunction-language derivation this leads to the 
extra term — Kr,dCr/2 in Eq. (I?7l) for Cr- This does not 
change the quasisteady value for Cf but changes the qua¬ 
sisteady value Cr = —igCe/[iArq + |t/|^/(xAfq + Kf/2) + 


Kr4/2], so that the Purcell rate is 


P = 2Re 


zA 


rq ~r |t/|^/(zAfq -|- K{l2) + /Ci._<j/2 


g^{Kq + Kr.d) 


A2 


(42) 

(43) 


instead of Eq. (15^ . Practically the same result can be 
obtained using the derivation via the density matrix evo¬ 
lution (assuming <C Kf). As expected, the dissipation 
Kr,d simply adds to the rate Kq seen by the qubit. Since 
Kr^d is not affected by the filter, it adds in the same way 
to the bandwidth Kr governing the qubit measurement 
process and thus deteriorates the Purcell rate suppres¬ 
sion dUD, replacing it with F = (Kq -I- Ki._d)/(Kr + K^^d)- 


IV. PURCELL RATE WITH MICROWAVE 
DRIVE AND BANDPASS FILTER 

The Purcell rate may decrease when the measurement 
microwave drive is added (^ . A simple physical reason 
is the ac Stark shift, which in the typical setup increases 
the absolute value of detuning between the qubit and 
readout resonator with increasing number of photons in 
the resonator, thus reducing the Purcell rate. However, 
this explanation may not necessarily work well quantita¬ 
tively. 

The Purcell rate suppression due to the microwave 
drive was analyzed in Ref. [2^ for the case without the 
Purcell filter and using the two-level approximation for 
the qubit. It was shown that in this case the suppression 
factor r(fi)/r(0) is approximately [(1 -I- h/ncrit)~^^^ + 
(1 -I- h/ucrit)”^]^ instead of the factor (1 -I- h/ucrit)”^ ex¬ 
pected from the ac Stark shift, where h is the mean num¬ 
ber of photons in the resonator and ricrit = (Arq/S^)^. 
This difference results in the ratio 3/2 between the cor¬ 
responding slopes of r(h) at small h, with the ac Stark 
shift model underestimating the Purcell rate suppres¬ 
sion (see the blue lines in Fig. 0]). However, when the 
third level of the qubit is taken into account, then the 
ac Stark shift model describes correctly the slope of r(h) 
at small h when the qubit anharmonicity is relatively 
small, [dq/Arql <C 1 (see Appendix). In this case the ac 
Stark shift model predicts r(h)/r(0) = 1-1- 4hx(0)/Ai.q 
at h <C u-crit, where x(0) is the value of x at h = 0; note 
that |x(0)| <C [g^/Arql and x(0) < 0 when Aj-q > 0. 

With the filter resonator, we also expect that the ac 
Stark shift model for the Purcell rate suppression should 
work reasonably well, so that from Eq. dSSl) we expect 

^_ 

^ [OJ,- Wq.eff(h)]2{[c.f - UqMnW + («f/2)n ’ 

(44) 

where ujq^esin) is the effective qubit frequency, 
Wq,eff(h) = Wq -I- 2x(0)h if wc neglect dependence of x 
on n and the “Lamb shift”. Therefore, in a typical situ¬ 
ation when |wf — ujr\ |wr — Wql and Kf <C [wr — Wql, we 
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FIG. 4: The Purcell relaxation rate r(h) with a microwave 
drive, normalized by the no-drive rate r(0), with the filter 
(red solid curve) and without the filter (blue solid curve), as 
functions of the mean number of photons h in the readout res¬ 
onator. The numerical simulations used the two-level approx¬ 
imation for the qubit, for which h/|Arq/x(0)| = h/4ncrit. The 
dashed lines show the values expected from the model based 
on the ac Stark shift: r(h)/r(0) = (1 + h/ricrit)”^ with the 
filter and (1 -I- h/ucrit)”^ without the filter. The parameters 
used in the simulations are given in the text. 


expect the suppression ratio 

r(fi) ^ Wr - Wq 

r(0) CUr CUq^eff(^) 


(45) 


To check the accuracy of this formula numerically, we 
need to add into the Hamiltonian the terms describ¬ 
ing the drive and higher levels in the qubit (see Ap¬ 
pendix). However, the resulting Hilbert space was too 
large for our numerical simulations, so we numerically 
calculated r(h) using only the two-level approximation 
for the qubit. Using the rotating frame based on the 
drive frequency Wd [i.e., Hq = Wd(o'+(T- -I-a’^a-I-6^6)], we 
then obtain the interaction Hamiltonian 


Vd = Ardo^o + Afd6^6 -b AqdCT+tT_ -b ^(aV- -b aa+) 

+ Qa^b + Q*ab^ + ErO^ + e*a, (46) 

where Am = wJ'-Wd, Afd = Wf-Wd, Aqd = Wq-Wd, and 
for simplicity we assumed that the drive Er is applied to 
the readout resonator. (For the rotating wave approxi¬ 
mation we also need to assume that lAml, |Afd|, |Aqd|, 
|g|, \0\, and Iffr/^rl are all small compared with Wd.) Note 
that in the two-level approximation [wr — coa pff(h)l^ = 
A^q -b Ag^n = A2q(l -b h/ricrit)- 

We have numerically solved the full master equation 
with the Hamiltonian ()46ll . including the decay K{ of 
the filter resonator. As the initial state we use the ex¬ 
cited state for the qubit and vacuum for the two res¬ 
onators, |'(/')in = |e00). The Purcell rate is extracted from 
the numerical solution of Pee{t) by fitting — ln[/5ee(i)] 
with a linear function in the long time limit (still re¬ 
quiring 1 — Pee(t) <C 1). In the simulations we pump 


the readout resonator with the frequency Wd = 
and control fi by choosing the corresponding value of 
fTp. The value of h is calculated numerically; it is close 
to what is expected from the solution of the classical 

|g\ 

field equations when the n-dependence of x ujI- 
is taken into account: in the two-level approximation 
X(n) = -gV[^rqyiTVh7^] and = uj^ + x(fi). 

Since changes with h, we change Wd accordingly. 

The red solid line in Fig.|3]shows the numerical results 
for the Purcell rate suppression factor F(h)/F(0) as a 
function ofh = normalized by |Arq/x(0)|. Note that 
in the two-level approximation (which we used in the sim¬ 
ulations) h/|Ai.q/x(0)| = fi/4ncrit- In the simulations we 
have used = 100 MHz, = 36 ns, = 0.71ns, 

I2tt = = 6.8 GHz, and Wq/STr = 6 GHz. The 

blue solid line shows the numerical suppression factor 
for the standard setup (without the filter resonator), 
also in the two-level approximation for the qubit. We 
see a larger suppression for the case with the filter, as 
expected from the ac Stark shift interpretation and the 
fact that F oc with the filter, while F oc in 

the standard setup. However, comparison with the pre¬ 
diction of the ac Stark shift model (dashed lines) does 
not show a quantitative agreement. There is about 10% 
discrepancy for the slope of r(h) between the red solid 
and red dashed lines in the case with the filter. This is 
a better agreement than for the case without the filter 
(blue lines). 

Note that the numerical calculations have been made 
using only the two-level model for the qubit. It is possible 
that the agreement between the ac Stark shift model and 
the numerical results is much better if three or more levels 
in the qubit are taken into account. This is an interesting 
question for further research. 

Also note that in experiments, increase of the drive 
power often leads to decrease of the qubit lifetime, in¬ 
stead of the increase, predicted by our analysis (both 
with and without the filter). The reason for this effect 
is not quite clear and may be related to various techni¬ 
cal issues. Therefore, either suppression or enhancement 
of the qubit relaxation with the drive power may be ob¬ 
served in actual experiments. 


V. CONCLUSION 

In this paper we have discussed the theory of the band¬ 
pass Purcell filter used in Refs. for measure¬ 

ment of superconducting qubits. An additional wide- 
bandwidth filter resonator (Fig. [2]) coupled to the read¬ 
out resonator easily passes the microwave field used for 
the qubit measurement, but strongly impedes the prop¬ 
agation of a photon at the qubit frequency, which is far 
outside of the filter bandwidth. A simple way to quanti¬ 
tatively describe the operation of the filter is by noticing 
that the effective decay rate Kes of the readout resonator 
[Eq. ([5])] depends on frequency. Therefore, the measure- 
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merit is governed by a relatively large value [Eq. (fTOl) ] , 
which permits a fast measurement, while the Purcell re¬ 
laxation is determined by a much smaller value Kq seen 
by the qubit [Eq. (ITT]) ] . The ratio of these effective decay 
rates of the readout resonator gives the suppression fac¬ 
tor for the qubit relaxation [Eq. (IT^ j. The result for the 
suppression factor is similar to the result obtained in Ref. 
[l7l| using circuit theory (with a few minor differences). 

We have first analyzed the operation of the Purcell 
filter quasiclassically, and then confirmed the results us¬ 
ing the quantum approach. In the quantum analysis we 
have used two approaches: based on decaying wavefunc- 
tion and density matrix evolutions. While the Purcell 
effect is traditionally described using density matrices, it 
is actually simpler to use the approach based on wave- 
functions. The results of our semiclassical and two quan¬ 
tum approaches are very close to each other; however, 
they are not identical because the approaches use slightly 
different approximations. A simple and most physically 
transparent result for the Purcell rate is given by Eq. 

dSl. 

The Purcell rate of the qubit decay is further sup¬ 
pressed when a microwave drive is applied for measure¬ 
ment. The effect is similar to what was discussed in Ref. 

for the standard setup without the filter, and can be 
crudely understood as being due to the ac Stark shift of 
the qubit frequency, which increases the resonator-qubit 
detuning A^q. Since the Purcell rate with the filter scales 
with Arq crudely as [Eq. dS^ ] instead of in the 
standard setup, the Purcell rate suppression due to mi¬ 
crowaves is stronger in the case with the filter. Numer¬ 
ical results for the suppression due to microwaves using 
the two-level approximation for the qubit show that the 
explanation based on ac Stark shift works well, but un¬ 
derestimates the effect by about 10%. This discrepancy 
may be significantly less if more levels in the qubit are 
taken into account; however, this still remains an open 
question. 

The bandpass Purcell filter decreases the qubit decay 
due to the Purcell effect by a factor of ^50 for typical 
parameters, for the same measurement conditions as in 
the standard setup without the filter. This allows much 
faster and more accurate measurement of superconduct¬ 
ing qubits, compared to the case without the filter. The 
qubit measurement within 140 ns with 99% fidelity using 
this filter has been demonstrated in Ref. [ 13 . With a 
slight change of parameters it seems possible to perform 
qubit readout within ^50 ns with fidelity approaching 
99.9%. Such fast and accurate qubit readout would be 
very useful for quantum information processing with su¬ 
perconducting qubits. 
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Appendix A: Qubit measurement and Purcell effect 
in the standard setup 

In this Appendix we review the trade-off between the 
Purcell rate and measurement time for a transmon or 
Xmon qubit in the standard cQED setup (without a fil¬ 
ter). In the standard setup [ii|,[l3[^ (Fig. H]) the qubit 
is dispersively coupled with the resonator, so that the 
qubit state slightly changes the resonator frequency. This 
change causes a phase shift (and in general an amplitude 
change) of a microwave field transmitted through or re¬ 
flected from the resonator. The transmitted or reflected 
microwave is then amplified and sent to a mixer, so that 
the phase shift (and amplitude change) can be discrim¬ 
inated, thus distinguishing the states \g) and |e) of the 
qubit. 


1. Basic theory 

For the basic analysis of measurement [2^, it is suf¬ 
ficient to consider three energy levels of the qubit: the 
ground state [g), the first excited state |e), and the sec¬ 
ond excited state [/), so that the Hamiltonian is (fi = 1) 

H = a;^|e)(e| -k {2uj^ - Jq)|/)(/| -k uj^a'^a 

+ ga'^\9){e\ + g*a\g){e\ + ga'^\e)(f\ + g*a\f)(e\ 

+ + H^ + (Al) 

where wj) is the bare qubit frequency, Jq is its ahnar- 
mononicity (Jq > 0, Jq/Wq 0.2 GHz/6 GHz <C 1), 
is the bare resonator frequency, a is the annihilation op¬ 
erator for the resonator, g is the coupling between the 
qubit and the resonator, g r; g is the similar cou¬ 
pling involving levels |e) and [/), and Sr is the normal¬ 
ized amplitude of the microwave drive with frequency Wd. 
For brevity describes the coupling of the resonator 
with the transmission line, which causes resonator en¬ 
ergy damping with the rate k, while describes the 
intrinsic qubit relaxation (excluding the Purcell effect) 
with the rate T’/’inj. Note that in the Hamiltonian (EH) 
we neglected the coupling terms creating or annihilating 
the double excitations in the qubit and the resonator. 
For simplicity we assume a real coupling: g* = g and 
9 * = 9 - 

For a simple analysis of the measurement process, let 
us first neglect Hj, and the drive Sr, and consider 
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the three Jaynes-Cummings ladders of states \g, n), |e, n), 
and I/, n), where n denotes the number of photons in 
the resonator. The coupling g provides the level repul¬ 
sion between \g,n + 1) and |e,n) (effective coupling is 
\/n+ 1 g), while g provides the level repulsion between 
\e,n+ 1) and |/, n) (with coupling + Ig). Assuming 
sufficiently large level separation, jwq— V^\g\ and 
Iwq — i5q — wj*! -y/n |g|, we can treat the level repulsion 
to lowest order; then the energies of the eigenstates \g, n) 
and |e, n) are 


\g,n) ■■ wb-a;b’ 

u, u in -\- 1) 

E] -r = w!? -b nw'? H-r-7— 

\e,n) q r b _ q.b 


ng^ 


— (5q — w. 


(A2) 

b- (A3) 


Therefore the effective resonator frequency when the 
qubit in the ground state is 


two-level case. This is because the transmon and Xmon 
qubits are only slightly different from a linear oscillator, 
for which g — ^/2g and 5q = 0, thus leading to x = 
0. The effect of nonzero anharmonicity (5q in Eq. (IA8D 
is more important than the effect of nonzero g — y/2g 
because 





(A12) 


and |A| <C Wq. Therefore x can be approximated as 


g^<5q ... g^^q 

A(A - <5q) A2 ’ 


(A13) 


where the last formula also assumes |A| (5q. 

The dispersive approximation (lAlOl) is based on the 
approximate formulas (IA2I) - (IA3|) for the eigenenergies, 
and therefore requires a limited number of photons n in 
the resonator, 


where 


= E 


|g,n-|-l) 


- El -r = 

Is.") 




(A4) 


A = Aq, = (A5) 

while for the qubit state |e) the effective resonator fre¬ 
quency is 


w, 


= Et 


^ - Wr + ^ - 


g 


A-5q- 


(A6) 


^|e,n-|-l) 

Denoting the frequency difference by 2x, we obtain 

= 2x, (A7) 


9 

A 


5V2 




g^-~gV‘^ 


A-^q A(A-^q) A-5q 

The corresponding effective qubit frequency is 


. .eff _ ZT- IT- . . .h \ y 

^--T — \ “I” 

q |e,n) \g,n) q ^ 


2nx, 


(AS) 


(A9) 


which includes the “Lamb shift” g^/A and the “ac Stark 
shift” 2nx- 

Therefore, in this case the first two lines of the Hamil¬ 
tonian (EH) can be approximated as 


n < min(ncrit, hcrit), '^crit 


A" . (A - 5q)2 

V’ 4^2 


(A14) 

where the factor 4 in the definitions of the critical photon 
numbers Ucrit and hcrit is a usual convention. For n be¬ 
yond this range it is still possible (at least to some extent) 
to use the dispersive approximation (lAlOl) if the spread 
of n is relatively small; however Wq, Wr, and x should be 
redefined. For that we need to use similar steps as in 
Eqs. (IA4I) - (IA9|) . but starting with more accurate formu¬ 
las than (IA2I) and (IA3|) . In particular, in the two-level 
approximation for the qubit (g = 0 oi Sq = oo) we would 
obtain x ~ 9^! \J^ + Ang'^ with h being the average 
(typical) number of photons; however, the generalization 
of the realistic case (IA13I) is not so simple (see below). 

The dispersive approximation (lAlOl) cannot reproduce 
the Purcell effect [llj after including the last line of 
the Hamiltonian (lAll) . so it should be added separately. 
Without the microwave drive = 0) only levels |e, 0), 
Ig, 1) and |g, 0) are involved into the evolution described 
by Eq. (ED, and for sufficiently small resonator band¬ 
width, K <C yA^ -I- 4g^, the qubit relaxation rate due to 
Purcell effect is 


r 


K 

2 



|A| \ 

a/A 2 -t4g2y 



(A15) 


H = 


UJo 


(Tz + uj^a)a + xa^a cr^. 


(AlO) 


where 


_ b , 5 
Wq — Wq -I- 


A’ 


LOr = Ulr — 


9 


2 A-5q 


(All) 


CTz = |e)(e| — |g)(g|, we shifted the energy by —a;q/2, and 
we no longer need the qubit state |/). 

Note that the dispersive coupling x given by Eq. (jA8l) 
is much smaller [ 23 | than the value g^/A expected in the 


Note that this rate does not depend on the qubit an¬ 
harmonicity 6q, in contrast to the dispersive coupling x, 
which vanishes at Jq —>■ 0. This is because the Pur¬ 
cell effect is essentially a linear effect (energy decay via 
decay of a coupled system), in contrast to x, which is 
based on qubit nonlinearity. This linearity is the rea¬ 
son why the Purcell rate P does not change (in the 
first approximation) when the microwave drive cre¬ 
ates a significant photon population in the resonator, 
1 n <C min(ncrit,hcrit). (Someone might naively ex¬ 
pect that P scales with n because of effective coupling 
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\/ng.) _ However, in the next approximation F depends the rate r|e)_>.|/) = K\{f,n — 2|a|e,n)p. These rates (up 

on n [23| and can be calculated as to sixth order in coupling) are 


T{n) = K\{g,n\a\e,n)f, (A16) 

with subsequent replacement of n with the average pho¬ 
ton number h if the spread of n is relatively small. 

Using the third-order (in g and/or g) eigenstates, 


ls)^k) 


Kg^n{n — 1) 




A2 


n 2 


Kg 

A2 


’^crit 


A2 A(2A - (5q) 
16(2A- 5q)2’ 


(A21) 

(A22) 


| 5 ,n) = 1 - 


ng^ 

2A2 

1 - 


\9,n) -t- 


\/n(n - l)gg 
A(2A - 5q) 


\f,n-2) 


Vngf^ ing"^ {n - l)g'^ 


2A2 A(2A-(5q) 


- 1 ), 


(A17) 


^ V 2A2 2(A-(5q)2y) 

VnTTg f 3(n-hl)5^ ng2(A - 2(5q) 

A ^ 2A2 2A(A - 5q)2 

x|g,n-H)- ^^^ |/,n-l), (AI 8 ) 


where the last term in Eq. (IA18|) for brevity is only of 
the first order, we find the Purcell rate 


r(ri) = ( 1 - % - 6n^ 


A2 


A2 


g2(3A - 4(5q) 
' A(A - 5q)2 


(A19) 


which is an approximation up to 5th order in g. This 
result coincides with the result of Ref. [ 2 ^ when g = 0. 
Approximating g = \/2 g [see Eq. (|A12p ] , we obtain 


r(n) 


2g^^q(2A- 3Jq) \ 

'"a2 V A2 A2(A-(5q)2 )■ 


(A20) 


It is interesting to note that while in the two-level ap¬ 
proximation (g = 0) the Purcell rate is always suppressed 
with increasing n [23| . Eq. (IA20p shows the suppression 
only when A < (3/2)(5q (which is the usual experimen¬ 
tal case, since A is usually negative). Numerical results 
using Eq. (IA16I) show that even when r(n) initially in¬ 
creases with n, it is still suppressed at larger n. Note 
that the result (|A20I) requires assumption (IA14I) of a 
sufficiently small nonlinearity. Comparing Eqs. (IA20I) 
and (IA13p . we see that in the case of large detuning, 
|A| (5q, the dependence of the Purcell rate on n is con¬ 

sistent with the explanation based on the ac Stark shift, 
r(n) ~ Kg2/(A -1- 2nx)^- (This is in contrast to the two- 
level approximation, in which this explanation leads to 
a discrepancy in the slope [ 2 ^ by a factor of 3/2.) We 
have checked that taking into account the fourth level in 
the qubit does not change Eqs. (IA19I) and (IA20I) : there 
are no additional contributions of the order g^. 

Besides the qubit energy relaxation F, the resonator 
damping k in the presence of drive leads to the qubit 
excitation [ 2 ^ |g) —>■ |e) with the rate F|g^^|g^ = 
K |(e, n — 2|a|g, n)p and the excitation |e) —>■ [/) with 


„ Kg’^nin — 1) f g^ 

^ (A - (5q)4 ^A-Jq “ 2A-5q 


2A-3<5q 

Kg2 (n/ncrit)^ ^qA® 

^ 2(A - (5q)6(2A - (5q)2(2A - 35q)2 ’ 


(A23) 

(A24) 


where g ~ V^g is the coupling due to the fourth qubit 
level with energy 3a;q — 3(5q, and in Eqs. (jA22l) and (IA24I) 

we also used g = V^g and replaced n{n — 1) with n^. 
In the assumed range, n <C ngrit, the excitation rates are 
much smaller than the Purcell decay rate F. 

Now let us discuss the n-dependence of the dispersive 
coupling x(n) = [ujf^n) — wi®^(n)]/2, where uji^\n) = 
E-, --rr — E -.—r and (n) = E-, --rr — E -,—r [see Eqs. 

|5,n+l) -ivy |e,n+l) |e,n) L ^ 

(ita . (iMi) . and dMI]. Using the three-level approxima¬ 
tion for the qubit with accuracy up to the fourth order 
in g and/or g (assuming n <C n-crit), we obtain 


(n) 


'a 


g‘^{2n + l) 2g'^g'^n 


9 


r_ 

A A-5„ 


A® A2(2A-5q)’ 

g‘^{2n + Z) g4(2n-fl) 


(A25) 


A® 


(A-<5q)3 


2Sqg'^g'^(n -h 1) 
A2(A-5q)2 


(A26) 


so that assuming g = V2g, we obtain 


X{n) 


g^^q 4g^l5q 1 - ^q -K ^q/2 

A(A-(5q)^ A4 (l-5q)3 


3ng4 1-52 + 53 _j 4/3 
^ A3 (i_^^/2)(l_5q)3’ 


(A27) 


where 6q = Sq/A. This result predicts a quite strong 
n-dependence of y, which is, however, not correct. The 
reason is that it is not sufficient to consider three qubit 
levels for x(n). When the fourth level of the qubit is 
taken into account (with energy Scu^ — 36q and coupling 

g), this does not change Eq. (|A25|1 for but intro¬ 

duces an additional term 


2g^g^n 

(A-5q)2(2A-35q) 


(A28) 
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into Eq. (IA26I) for u}i‘^\n). Assuming g = VSg, this 
changes Eq. (IA27I) into 


Xin) Ki - 


9% 


A(A - (5q) 


4g^l5q 1 - ^q + ^q/2 

A" (1 - 5q)3 


9ngX 1 - (5/3)Jq + (11/9)J^ - 
2A5 (1-V2)(l-Jq)3(l-3V2)' 

which shows a quite weak dependence on n, 


(A29) 


dxjn) _ 9 ^ x(0) 

dn 8 A ricrit ’ 


(A30) 


assuming Jq ^ |A|. In the usual case when A < 0, the 
absolute value of x(’^) decreases with increasing n. 

Note that in Eqs. (IA27I) and (IA29I) we used g = \/2g 
and g = \/ig. If a better approximation is used, 


g Ri y2g(l - (5q/2a;q), g Ri VSg (I - Jq/Wq), (A31) 


then correction to the g^-term in Eq. (IA26I) creates an 
additional contribution g^(5q/[u;q(A — Jq)] to X: which for 
typical parameters is much larger than the second terms 
in Eqs. (IA27I) and (IA29|) . The correction (IA3II) also yields 
an additional n-dependent contribution of approximately 
2g^5qn/(A3a;q) (assuming 5q <C |A|) to both ujl^\n) and 

uji^\n). These additional slopes are comparable to the 
slope of xin) in Eq. (IA29|) because A^/((5qa;q) is on the 
order of 1 for typical experimental parameters. However, 
the contribution to x(^) from the correction (IA3II) is only 
9g'^5qn/(A^a;q) (assuming Jq <C |A|), which is smaller 
than the last term in Eq. (IA29|) because |A| <C Wq. 
Therefore, the slope of xin) in Eq. (IA29I) is practically 
not affected by the correction (IA3II) . Similarly, inaccu¬ 
racy of our approximation of the fourth qubit level energy 
by dWq — 3(5q produces only a small correction to the slope 
of xin)- the correction to the fourth-level energy is on the 
order of ^q/w^, and therefore the correction to xin) [via 
Eq. llA28|) ] is on the order of g"‘^qn/(A'*Wq) for 5q <C |A|. 

Thus, we see that for a transmon or Xmon qubit, cal¬ 
culation of x(0) requires at least three qubit levels to 
be taken into account, while the first correction due to 
xin) dependence (n <C Ucrit) requires at least four qubit 
levels to be considered. In contrast, calculation of the 
Purcell rate r(0) requires two qubit levels, while the first 
correction due to r(n) dependence requires three qubit 
levels. 


2. Measurement error 

Now let us discuss the qubit measurement error caused 
by the qubit relaxation due to Purcell effect. To describe 
measurement, we will use the dispersive approximation 
(IA10|) with X = —g^dq/A^ [Eq. (IAI3I) ] and neglect the 
small difference between A defined in Eq. (IA5I) and Wq — 
Wr defined via Eq. (lAllI) . For the Purcell rate we will 


use the simplest form P = g^re/A^. Both approximations 
assume a sufficiently small number of photons, Eq. (IAI4I) . 

Assuming that the qubit is either in the state |g) or 
|e) (non-evolving), from Eq. (lAIOl) we see that the ef¬ 
fective resonator frequency is constant in time, Wr ± Xj 
where the upper sign is for the state je) and the lower 
sign is for jg). The corresponding resonator state is 
a coherent state |Qf±), characterized by an amplitude 
a±it) in the rotating frame based on the drive fre¬ 
quency Wd [in the lab frame the resonator wave func¬ 
tion is overall 

phase]. The evolution of the amplitude is 

d± =-f(Ard ± x) a± -Q!± - lEr, (A32) 

where the upper sign everywhere is for the state | e), 

Ard = Wr-a;d, (A33) 

and complex Srit) can in general depend on time to de¬ 
scribe a drive with changing amplitude and frequency. 
(Note a difference in notation compared with the main 
text: now A^d does not depend on the qubit state and 
the frequency shift ±x is added explicitly.) The steady 
state for a steady drive, Sr = const, is 


a± 


—ie-c 

k /2 + f(Ard ± x) ’ 


(A34) 


so that the two coherent states are separated by 


_ -2erX 

(«/2 + *Am)"+x'’ 


(A35) 


and this is the difference, which can be sensed by the 
homodyne detection (it does not matter whether the mi¬ 
crowave is transmitted or reflected from the resonator 
when we discuss the measurement in terms of a±). 

In the case when A^d = 0 (so that Wd = Wr), both 
states have the same average number of photons, n = 
|a+P = jo-P = |£r|^/(K^/4-|-x^), and the absolute value 
of the state separation is 


(5Q = |a+-a-|= 

\/(k/2x)^ + 1 

Theoretically, the distinguishability of two coherent 
states separated by 5a is the same as the distinguisha¬ 
bility of two Gaussians with width (standard deviation) 
of 1/2 each, separated by 5a. However, a measurement 
for the duration increases the effective separation by 
a factor of \/ntm, while imperfect quantum efficiency t] 
of the measurement (mainly due to the amplifier noise) 
increases the Gaussian width by or, equivalently, 

decreases the effective separation by Therefore, 

we may think about the distinguishability of two Gaus¬ 
sians with width 1/2 each, separated by 

^^eft — \/g^^m 5a^ 


(A37) 
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which gives the error probability due to state separation 

psep ^ 1 - Erf((5aeff/y2) _ exp[-(foeff)V2] 

"" 2 V^foeff ■ ^ 

The other contribution to the measurement error Perr 
for the state |e) comes from the qubit energy relaxation 
with the rate T TTfint during the measurement time tm, 

Perr = P^rT + ^ + T’urrt), (A39) 

where the factor 1/2 is because the relaxation moment 
is distributed practically uniformly within the measure¬ 
ment duration tm- Since P|/p decreases with time tm 
exponentially, the main limitation for Perr comes from 
the second term. 

It is easy to find that P|/p = 10“^ corresponds to 
5aefi = 2.3, P|/P = 10“^ corresponds to (Joeff = 3.1, 
and P|/r^ = 10“'^ corresponds to Jaeff = 3.7. For an 
estimate let us choose Saes > 3. Then from Eq. (IA37I) 
tm > 9/{r]KSa^), and therefore even neglecting intrinsic 
relaxation r/"jnt in Eq- (IA39I) . we obtain the condition 

E < iperr?7«(<^a)^- (A40) 

Now using r = Kg'^jA^ and assuming k > 2 Ixl in Eq. 
(IA36I) . so that 6a ~ AxV^Ik, we rewrite this condition 
as 


Ixl 


l5^q 


K < 2v^Perr?7« ^ - 2y/PerrVn , 


(A41) 


where for the second expression we used x — 

Since the measurement time t„i should be at least few 
times longer than k~^, we obtain 


4 ^ 2 |A| 

f > _ > _ . J_L 

“ K y/PerrW IS'^ql 


4/^q 

\/ Perr^'x/ u/Ucrit 


(A42) 


These estimates show that the Purcell effect requires 
a sufficiently long measurement time when we desire a 
small measurement error Perr- The limitation is not se¬ 
vere, but it is still inconsistent with a fast accurate mea¬ 
surement needed for quantum error correction. As an ex¬ 
ample, from Eq. (|A42ll we see that for 5q/27r ~ 200 MHz, 
Perr — 10“^, T] ~ 0.3 and h ~ ncrit/4 we need > 400 ns 
(this in turn would require a very long Ti^nt)- Also, the 
detuning should be sufficiently large. 


\A/g\ > v'Ktm/2Perr > (A43) 


as directly follows from T = /c(g/A)^ and Kfm > 4. 

As a more detailed example, let us choose g/2Tr ~ 
30 MHz, A/27r ~ -1.35 GHz, SJ2tt ~ 200 MHz, k"! ~ 
100 ns, tin — 400 ns, g ~ 0.3, and n ~ ncrit/4 — 125. 
Then for the excited qubit state the Purcell decay brings 
the error contribution tmF/2 ~ 10“^; the dispersive cou¬ 
pling is x/27r ~ —0.1 MHz, so Sa ~ 0.25-\/fi ~ 2.8 and 
(5aeff — 3, so the separation error is also about 10“^ [actu¬ 
ally slightly larger since |x| decreases with n for negative 
A - see Eq. (IA30|) ]. 

Thus, we see that even with the qubit decay due to 
the Purcell effect, it is possible to measure a qubit with 
a low measurement error, but this requires a relatively 
long time and large resonator-qubit detuning. Note that 
we also had to assume a large number of photons in 
the resonator, which can lead to detrimental effects (ne- 
glected in our analysis) such as dressed dephasing [32Ll33l| 
and various imperfections related to nonlinear dynam¬ 
ics. Suppression of the Purcell effect (using a Purcell 
filter or other means) allows us to significantly increase 
the ratio [^/Aj, thus increasing the coupling x and cor¬ 
respondingly decreasing the measurement time for the 
same measurement error. 
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